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Abstract 

In this paper we study the long time asymptotic behavior for a class of diffusion- 
aggregation equations. Most results except the ones in Section 3.3 concern radial so- 
lutions. The main tools used in the paper are maximum-principle type arguments on 
mass concentration of solutions, as well as energy method. For the Patlak-Keller-Segel 
problem with critical power m = 2 — 2/d, we prove that all radial solutions with critical 
mass would converge to a family of stationary solutions, while all radial solutions with 
subcritical mass converge to a self-similar dissipating solution algebraically fast. For non- 
radial solutions, we obtain convergence towards the self-similar dissipating solution when 
the mass is sufficiently small. We also apply the mass comparison method to another 
aggregation model with repulsive-attractive interaction, and prove that radial solutions 
converge to the stationary solution exponentially fast. 



1 Introduction 

Recently there has been a growing interest in the study of nonlocal aggregation phenomena. 
The most widely studied models are the Patlak-Keller-Segel (PKS) models, which describes 
the cell movement driven by chemotaxis jKSilPllH] . In this paper, we study the PKS equation 
in dimension d > 3 with degenerate diffusion, given by 

Ut = Aw" -t- V • (uVc) X G M'', t > 0, 

< Ac =u X eM.'^,t>0, (1.1) 

^ u{x,0) =ua{x) e Ll{R'^:{l + \x\^)dx)nL'^{m.'^) xeM.'^, 

here m > 1 models the local repulsion of cells with anti-crowding effects |TB1 ITBL] . This 
model is a generalization of the classical parabolic-elliptic PKS model in 2D, which has been 
extensively studied over the years (see the review [Hl|B] and |DP[ IBDPl IBCMj ). 

Throughout this paper, we focus on the critical power m — nid 2 — 2/(i, which produces 
an exact balance between the diffusion term and the aggregation term when one performs 
a mass- invariant scaling. To study the wcU-poscdness of (|l.ip . the following free energy 
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functional (jl.2[) is an important quantity, where the first term is usually referred to as the 
entropy and the latter term is referred to as the interaction energy: 

T[u) = I ( —^vJ^ + -u{u * M)] dx. (1.2) 

Jri \m - 1 2 J 

Here A/'(x) = " {d-\)a-a l-'^P"'^ Newtonian potential in d > 3, with (7^1 the surface area of 

the sphere S''"^ in K**. The free energy of a weak solution to (jl.ip is non-increasing in time; 
indeed, it is shown that (jl.ip is the gradient flow for T with respect to the Wasserstein metric 
(see for example [AGS] and fCMV] ). 

The key observation in jBCL] is the sharp Hardy-Littlcwood-Sobolev inequality, which 
bounds the interaction term in the free energy functional by the entropy term: 



u{u * N)dx 



<C^*hll?/wJI«lir"m.v (1-3) 



where C* is a constant only depending on the dimension d. Making use of this inequality, it 
is proved in |BCLj (and generalized by |BRB| for more general kernels) that there exists a 
critical mass only depending on d, which sharply divides the possibility of finite time blow 
up and global existence. Although the global existence/blow-up results arc well understood, 
the asymptotic behavior of the solution has not been fully investigated: this motivates our 
study. In this paper we study the asymptotic behavior of solutions with mass M < M^, and 
try to answer the open questions raised in |BCL| and [B] . 

o Critical case: 

When M = Mc, it is proved in |BCLj that the global minimizers of the free energy functional 
J-' have zero free energy, and are given by the one-parameter family 

1 X 

ur{x) - -^M^) (1-4) 

subject to translations. Here ui is the imique radial classical solution to 

^^Au'"-^ +ui = in B(0, 1), with ui = on dB(0, 1). (1.5) 
m — 1 

It was unknown that whether this family of stationary solutions attract some solutions. 

In Section 3.1, we prove the convergence of radial solutions towards this family of stationary 
solution (see Theorem 13. 2p using a combination of mass comparison and energy method. The 
mass comparison property is a version of comparison principle on the mass distribution of 
solutions (see Proposition 12. 4p . which has been introduced in jKYj . Although it only works 
for radial solutions, it provides more delicate control than the energy method: it has been 
recently used in |BK| to prove finite time blow-up of solutions with supercritical mass for a 
diffusion-aggregation model where Virial identity does not apply. We mention that the mass 
comparison property have been previously observed for the porous medium equation jV] and 
PKS models with hnear diffusion [JLllBKLNj . 

o Subcritical case: 

When < Af < Mc, the weak solution exists globally in time, as long as its initial L"'-norm 
is finite |BCL1 IBRBj . Moreover it has been proved in jBCL| that there exists a dissipating 
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self-similar solution, with the same scaling as the porous medium equation. However it was 
unknown whether this self-similar solution would attract all solutions in the intermediate 
asymptotics. 

In Section 3.2, we prove algebraic convergence of radial solutions towards this self-similar 
dissipating solution, given that the initial data is bounded and compactly supported. This is 
done by constructing explicit barriers in the mass comparison sense. 

For general (non-radial) solutions, the asymptotic behavior of solution to is unknown 
for all mass sizes when m = 2 — 2/d. In Section 3.3, we prove that when the mass is 
sufficiently small, every solution to with compactly supported initial data converges to 
the self-similar dissipating solution. We mention that similar results for small-mass solutions 
are obtained by |BEDFj for the PKS equation in 2D with linear diffusion, however their 
argument is based on a spectral gap method and cannot be generalized to (jl.ll) due to the 
nonlinear diffusion term. 

When dealing with non-radial solutions with small mass, our key result is the uniqueness of 
stationary solution (in rescaled variables), which is proved using a maximum principle type 
argument (see Theorem 13. lip . We point out that although it is well known that the global 
minimizer to the rescaled energy p.l5[) is unique |BCL| . there are few results concerning the 
uniqueness of stationary solutions for nonlocal PDEs, except in the following special cases: 
the stationary solution to a similar equation is proved to be unique by |BDF| in the ID case, 
and stationary solution to the 2D Navier-Stokes equation (in rescaled variable) is proved to 
be unique by |GGL| , where their proof is based on the fact that all stationary solutions have 
the same second moment, thus cannot be applied to (jl.ll) . 



In Section 4 we generalize the mass comparison methods to an aggregation equation 

ut - V ■ (wVif * u) = 0, (1.6) 

and prove that the radial solution converges towards the stationary solution exponentially 
fast for a class of kernel K. Equation (|1.6p appears in various contexts as a mathematical 
model for biological aggregations |ME1 ITBLj . In order to capture the biologically relevant 
features of solutions, it is desirable for the interaction kernel K to be repulsive for short-range 
interactions and attractive for long-range. Aggregation equations with such kinds of kernel 
are studied in jBCLRl HTBI IFHKl iFHl IKSUB] . In this paper we adopt the kernel K proposed 
in |FHK| ■ which has a repulsion component in the form of the Newtonian potential M and 
an attraction component satisfying the power law: 

K{x)^M{x) + -^\x\\ (1.7) 

here when q — Q the second term is replaced by ln|x|. We assume that the attraction part 
is less singular than the repulsion part at the origin, i.e. q > 2 — d. In addition, we assume 
q < 2, i.e. the long-range attraction does not grow more than linearly as the distance goes to 
infinity. Note that in |FHj they mostly focus on the case q >2. 

The existence and uniqueness of weak solution is established in jFHj for q > 2 — d. In 
addition, they proved that for any mass there is a unique radial stationary solution which is 
continuous in its support. While numerical results in |FH| suggests that this stationary solu- 
tion should be a global attractor, there is no rigorous proof. We prove that when 2 — d< q < 2, 
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the mass comparison argument can be easily generalized to (jl.6|) . then we construct barriers 
in the mass comparison sense to prove that all radial solution converges to the stationary 
solution exponentially fast. 

Outline of the Paper. In Section 2 we state a mass comparison result for a general diffusion- 
aggregation equation. In section 3 we apply it to the PKS model with critical power, and 
obtain some asymptotic results for radial solutions with critical and subcritical mass sizes 
respectively. In Section 4 we demonstrate that the mass comparison can also be applied to an 
aggregation model with repulsive-attractive interaction, and prove the asymptotic convergence 
of solution towards the stationary solution. 

1.1 Summary of Results 

By constructing explicit barriers in the mass comparison sense and using energy method, we 
obtain the following results for radial solutions of 

Theorem 1.1. Suppose d > 3 and m = 2 — 2/d. Let u{x,t) he the weak solution to (jl.ip 
with mass A and initial data uq G L^{M.'^; (1 -|- H L^iW^), where uq is continuous, 

radially symmetric and compactly supported. Then the following results hold: 

• If < A < Mc, u{-,t) converges to the dissipating self-similar solution ua as t ^ oo, 
where the Wasserstein distance between t) and ua decays algebraically fast ast oo. 
(Corollary\3lM> 

• If A = Mc and uq satisfies VuJJ^ € L^(R'') in addition to the assumptions above, then 
u{-,t) ujig in L°°{M.'^) as t ~^ oo for some Rq > 0, where uji^ is a stationary solution 
defined in (TheoremW^) 

For general (possibly non-radial) initial data, we use a maximum principle type method to 
prove that when the mass is sufRciently small, every compactly supported stationary solution 
must be radially symmetric. This leads to the following asymptotic convergence result: 

Theorem 1.2. Suppose d > 3 and m — 2 — 2/d. Let u[x, t) be the weak solution to with 
mass < A < A/c/2 being sufficiently small, and the initial data uq G L^(R'^: (1 + |a:p)da;) H 
L°°(R'*) is continuous and compactly supported. Then we have 

lim \\u{-, t) — ua\\p = for all I < p < oo, 
where ua is the self-similar dissipating solution defined in p.l7p . f C'orollarv \3.LS\) 

In Section 4 we generalize the mass comparison methods to the repulsive-attractive aggrega- 
tion equation (|4.ip . and obtain the following asymptotic convergence result ior 2 — d < q < 2: 

Theorem 1.3. Assume 2 — d < q < 2. Let u be a weak solution to (|1.6p with initial data uq 
and mass A, where Uq G L"'^(M'^) n L°°(M'^) is non-negative, radially symmetric and compactly 
supported. In addition, we assume that Uq is strictly positive in a neighborhood ofO. Let Ug 
be the unique stationary solution with mass A, as given by Proposition \4-. 2\ Then as t ^ oo, 
u(-,t) converges to Us exponentially fast in Wasserstein distance. (Theorem \4.5\ l 
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2 Mass Comparison for Radial Solutions 

In this section we consider the following type of diffusion-aggregation equation 

ut = ciAu™ + V • (uV(u * (caAA + C3/C) + V)), (2.1) 

where Af{x) = — (■^_2)g^ is the Newtonian potential in d > 3, with ad the surface area 

of the sphere S''"^ in Mf^. We make the following assumptions on the kernel /C, the potential 
V and the coefhcicnts: 

(C) ci, C3 > 0, and C2 G M can be of any sign. 
(Kl) /C is radially symmetric. 

(K2) A/C e Li(R^), A/C > and is radially decreasing. 
(VI) V e C^{R'^) is radially symmetric. 

For a radially symmetric function it(a;, t), we define its mass function M(r^ t; u) by 

M{r,t;u) := / u{x,t)dt, (2.2) 

JB{0,r) 

and we may write it as Al{r,t) if the dependence on the function u is clear. The following 
lemma describes the PDE satisfied by the mass function. 

Lemma 2.1 (Evolution of Mass Function). Let u[x,t) he a non-negative smooth radially 
symmetric solution to (j2.ip . Let M{r,t) = M{r,t;u) be as defined in (j2.2p . Then M{r,t) 
satisfies 

-—-^Ci(7dr dr{ -f-r) +drM [^drMdrV, (2.3) 

where M{r, t\ u) := /b(-q r) " * AlCdx. 

Proof. Due to divergence theorem and radial symmetry of u, wc have 

^ = ad/-l[cia,?/™ + U{dr{u * {C2N + C3/C)) + drV)]. (2.4) 



Note that radial symmetry of u also gives 



u{r) = (2.5) 
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It remains to write dr{u*J\f) and dr{u*lC) in terms of AI. For dr{u*J\f), divergence theorem 
gives 

^, lBiO,r)^^*J^dx M{r,t) 

ar(u*J\) = ^— ; = ^— r- (2.6) 

We can similarly obtain 

Or(u*/C)^ ^— ; = r-Ti 2.7) 

where M{r, t; u) := /^^g u * A/Cdx. Plug (|23|) . ((2^ and (fSj]) into equation ([2^ . and then 
we can obtain (12.31). □ 



Definition 2.2. Lef ui anrf U2 be two non-negative radially symmetric functions in L^(K''). 
We say ui is less concentrated than U2, or ui < U2, if 



ui{x)dx < / U2{x)dx for all r > 0. 

B(0,r) ■JB{0,r) 

Definition 2.3. Let ui{x,t) be a non-negative, radially symmetric function in L^{M.'^) 
L°°{M.'^), which is in its positive set. We say ui is a supcrsolution of p.ip in the mass 
comparison sense if Mi{r,t) := M{r,t;ui) is a supcrsolution of p.3p . i.e. Mi{r,t) and 
Mi{r,t) := M{r,t;ui) satisfy 

^rr ^ ciCTrfr" — -)™ + drMi — I + drMidrV, 2.8) 

in the positive set of Ui. 

Similarly we can define a subsolution of (|2.1[) in the mass comparison sense. 

Proposition 2.4 (mass comparison). Suppose in > 1, and ci, 02,03, V,IC satisfy the as- 
sumptions (C),(K1),(K2),(V1). Letui{x,t) be a supcrsolution and U2{x,t) be a subsolution 
of (|2.ip in the mass comparison sense for t G [0,r]. Further assume that both Ui are bounded, 
and Ui 's preserve their mass over time, i.e., J Ui{-,t)dx and J U2{-,t)dx stay constant for all 
< t < T. Then their mass functions are ordered for all times: i.e., if ui{x,0) >- U2{x,0), 
then we have ui{x,t) >- U2{x,t) for all t G [0,T]. 

Proof. Let Mi{r,t) be the mass function for Ui, where i = 1,2. We claim that Mi{r,t) > 
M2{r,t) for all r > and t £ [0,T], which proves the proposition. 

For the boundary conditions of Mi, note that 
fA/i(0,t) = M2{Q,t) = for ah t G [0,T], 

\\imr^^{Mi{r,t) - M2{r,t)) ^ Jj^,{ui{x,0) - U2{x,0))dx > for ah t e [0,r]. 
As for initial data, we have Mi{r, 0) > M2{r, 0) for all r > 0. 
For given A > 0, we define 

wir,t) := {M2ir,t)-Mi{r,t))e-^', 
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where A is a large constant to be determined later. Suppose the claim is false, then w attains 
a positive maximum at some point (ri, ti) in the domain (0, oo) x (0, T]. Moreover, since the 
mass of both ui and U2 are preserved over time and thus are ordered, we know that {ri,ti) 
must lie inside the positive set for both ui and U2, where Mi's are C'^'l- 

Since w attains a maximum at (ri,ti), the following inequalities hold at (ri,ii): 

> =^ dt{M2 - Ml) > \{M2 - Ml) (2.9) 

= =^ drMi = drAh > (2.10) 

Wrr <0 =^ drrMi > drrM2 (2.11) 



Now we will analyze the terms on the right hand side of (|2.8|) one by one. For the first term, 
(|2J0)) and ((2TTt imply that 

Ci5.(^4^)"-cia,.(A^)-<0 at(ri,ti). (2.12) 

For the term coming from Newtonian potential, we have 

^^^^^C2(M2_Mi) ^ _ < u^^M{M2 - Ml), (2.13) 

where Wmax ■= nia'X{supRdx[o,T] ""Ij supRdx[o,T] ""2} is finite by assumption on ui and U2- 
We next claim 

Q^Mi^^^^^^^^^^^^^^^ < u^,^C3\\AICUM2 - Mi){n,h). (2.14) 
To prove the claim, note that M2 — Mi can be rewritten as 

Ivhiruh) - Miin^ti) = / {{u2 - ui) * AlC) XB{o,r,)dx 

(U2 - Ml) (XB(0,ri) * AlCjdx. (2.15) 

Note that A/C > is radially decreasing due to assumption (K2), thus Xs(o,ri) * A.K. is non- 
negative, radially decreasing and has maximum less than or equal to ||A/C||i. Therefore we 
can use a sum of bump function to approximate Xs(o,ri) * A/C. where the sum of the height 
is less than ||A/C||i. Hence 

M2(ri,ti) -Mi(ri,ti) < || A/C||i sup(A/2 - Afi)(a;, <i) = || A/C||i(il/2 - Afi)(ri, ii), 

X 

which proves the claim (|2.14p . Finally, for the last term in (|2.8p , as a result of (|2.10p we have 

drM2dV - drMidV = 0. (2.16) 

Now we subtract (|2.8p with the corresponding equation for the subsolution. Due to the 
inequalities ([2TT2)) . ((2?T3| . ([2TT4)) and ((2ll| . we obtain that 

dt{M2 - Ml) < «„,ax(|c2| + C3||AAli)(M2 - Mi). 

Hence if we choose A > Umax(|c2| + C3||Aif||i) in the beginning of the proof, the inequality 
above will contradict (12.91). □ 
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Remark 2.5. // both ui and U2 are supported in some compact set B(0, R) for < t < T , 
then the assumption (K2) can he replaced by (K2') as follows: 

(K2') AX e L]^^{M.'^), AK > and is radially decreasing. 

The proof under condition (K2') is almost the same, except that there is some change in 
(|2.14p . Note that in this case (j2.15p still holds, and we can bound the maximum ofxB{o.ri)*'^^ 
by Jg^Q ^-j AlCdx. This yields an inequality similar to (|2.14p . with the ||A/C||i in the right hand 
side replaced by j^^^ AlCdx. 

3 Application to Critical Patlak-Keller-Segel Models 

3.1 Convergence towards stationary solution for critical mass 

In this subsection, we prove that every radial solution with mass Mc and continuous, com- 
pactly supported initial data will be eventually attracted to some stationary solution within 
the family pTi]) . 

If the initial data is bounded above and has the critical mass Mc, then it is proved in 
[BCLj that the weak solution to exists globally in time. In the next lemma we prove the 
solution has a global (in time) L°° bound. In addition, if the initial data is radially symmetric 
and compactly supported, then the support of the solution would stay uniformly bounded in 
time. 

Lemma 3.1. Suppose d > 3 and m — 1 — 2/ d. Consider the problem (jl.ip with initial data 
Uq € (1 + |a;p)d2;) D L'^{W^), where Uq is continuous and has critical mass M^- Then 

the L°° norm of the weak solution u{x, t) is globally bounded, i.e. there exists K > Q depending 
on ||mo||oo o,nd d, such that \\u{-, i)||ioo(j{d-) < Ki for all t > 0. 

// Uq is radially symmetric and compactly supported in addition to the assumptions above, 
then there exists some R2 > 0, such that {u{-,t) > 0} C i3(0,i?2) for all t >Q, where R2 
depend on d and uq. 

Proof. In order to bound the L°° norm of u(-,t), we first consider equation (jl.ip with 
symmtetrized initial data, which is described below. Let u{-,t) be the solution to with 
initial data Uq, where Uq{-) is the radial decreasing rearrangement of uq. Here the radial 
decreasing rearrangement of a non-negative function / is defined as 

f*{x):^ X{f>tY{x)dt. (3.1) 
Jo 

Since u has a radially symmetric initial data and has mass Mc, due to jBCLj . we readily 
obtain that u exists globally in time, and u is radially symmetric for alH > 0. We first prove 
that there is a global L°° bound for u. 

Since ||u(-,0)||oo = |[ito||oo < 00, we can choose i?i > depending on ||mo||oo, where Ri 
is sufHciently small such that < u ^ ^ , where u ^ ^ is as defined in (jl.4p . Then the mass 
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comparison result in Proposition 12.41 yields that 

u{-, t) -< UR^ for all t > 0. (3.2) 

Now we go back to the original solution u, and compare u with u. It is proved in Theorem 
6.3 of 112 that 

u*{-,t) -< u{-,t) for all t > 0. 
Combining the above two inequalities together, we readily obtain that 

^ Ufli for all i > 0, 

which implies that u*{Q,t) < ur-^{0) for t>0. Note that u*{-,t) is radially decreasing for all 
f > by definition, and ur^ is radially decreasing due to jBCLj . Hence the above inequality 
implies that 

\H-,t)\\^ = \\u*i-,t)\\^ < \\UR,\U = i?rV(0) for t > 0, (3.3) 
thus u has a global L°° bound i?j^'*ui(0), where ui is as defined in (|1.4I) . 

Next we hope to show that if uo is radially symmetric and compactly supported in addition 
to the conditions above, the support of u{-,t) will stay in some compact set for all time. 
We first prove it for the case where Mo(0) > 0. Due to the continuity of uq, we have Uq 
is uniformly positive in a neighborhood of 0. This enables us to choose i?2 > sufficiently 
large such that uq >- ur^ , where ur^ is as defined in (|1.4p . Proposition 12.41 again gives us 
u{-,t) >- ur2 for alH > 0, which implies that 

supp u{-, t) C supp UR^ = B{0, R2) for all t >0. 

If uo(0) = 0, we claim that after some finite time ti, u(0,ii) becomes positive, and u{-,ti) 
has a compact support, where ti depends on d and uq. Then we can take ti as the starting 
time and argue as in the case mo(0) > 0. 

Now we will prove the claim. This is done by performing the mass comparison between u 
and w, where w is the solution to the porous medium equation 

= Aw™ inK'* X [0,00), (3.4) 

with initial data w{-,0) = uq. It can be readily checked that m is a supersolution of (|3.4p in 
the mass comparison sense, hence Proposition 12.41 imphes that u{-,t) >- w{-,t) for all t > 0, 
which yields 

u{0,t) > w{0,t) and supp t) C supp w(-, t) for alH > 0. (3-5) 

For the porous medium equation (|3.4p . it is well known that the solution will eventually 
converge to the self-similar Barenblatt profile (see [V] for example), which has a positive 
density at for all t > 0. Hence there exists some ti > such that w{0,ti) > 0. In 
addition, w{-,ti) has a compact support, due to the finite speed of propagation property of 
porous medium equation jVj. Therefore (|3.5p yields that u{0,ti) > w{0,ti) > and u{-,ti) 
is compactly supported, which prove the claim. □ 



9 



Next we prove that under the conditions in Lemma 13.11 every radial sohition converges to 
some stationary solution in the family (jl.4[) as t — > oo. To do this we investigate the free 
energy functional (|1.2p , and make use of the following result proved in |BCLj and |BRB| : Let 
u be a weak solution to (|l.ip . then it satisfies the following energy dissipation inequality for 
almost every t during its existence: 



1 



V7V* ' 



dxdt < J-(uo). 



(3.6) 



Theorem 3.2. Suppose d > 3 and m = 2 
critical mass Mc and initial data uq G L]^ (. 



2/d. Let u{x,t) be the weak solution to (|l.ip with 
{l + \x\'^)dx)nL°° {R'^), where uq is continuous, 
radially symmetric and compactly supported, and satisfies Vu™ € L'^i'M.'^). Then there exists 



Rq > depending on uq and d, such that u(-,t) 
as defined in (11.41) . 



UR^, in L°°(M'^) as t 



where is 



Proof. Due to Lemma |3. 11 we obtain the existence of a weak solution globally in time, which 
has a global L°° bound. In addition, by treating u * A/" as an a priori potential in (jl.ip and 
applying the continuity result in [D], we obtain that u{x,t) is uniformly continuous in space 
and time in [t, c») for all r > 0. 

Our preliminary goal is to find a time sequence {tn} which increases to infinity, such that 
u(tn) uniformly converges to some stationary solution as n — > oo. Note that J-{u{-,t)) is 
non-increasing for almost every t due to (j3.6p . and is bounded below as i — > oo. This enables 
us to find a time sequence {tn} which increases to infinity, such that 

2 



lim / u{tn) 



m — 1 



Vu(t„)'"-i + VA/'*u(t„) 



dx = 0. (3.7) 



We will slightly abuse the notation and denote u(tn) by u„. Note that is uniformly 

bounded and equicontinuous, hence Arzela-Ascoli theorem enables us to find a subsequence 
of {un}, such that 

Un ^ ^oo 

uniformly in n, (3.8) 

where Uqo is some radially symmetric and continuous function. Moreover, Lemma l3 . 1 1 implies 
that the support of all stays in some fixed compact set, hence we have Uoo is compactly 
supported as well, and it has mass M^- We will prove that Uoo is indeed a stationary solution 
later. 

Wc next claim that {Vu™} arc uniformly bounded in L'^{W^). To prove the claim, note 
that 



|Vm™ + u„\7J\f * Unl'^dx < 



2 



dx. 



where the right hand side is uniformly bounded for all n. In addition, since are uni- 

formly bounded and are all supported in some B{0, R), wc know J^^ m„| VA/" *Un\'^dx is also 
uniformly bounded for all n. Therefore triangle inequality yields the uniform boundedness of 
/jjji |VM™p(ix, which proves the claim. 

The uniform boundedness of {Vu™} in i^(R'^) implies that {Vit™} is in i^(]R'') as well. 
Moreover, we can find a subsequence of {u„} (which we again denote by {u„} for the simplicity 
of notation), such that 
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V<' ^ Vu'^ as n ^ oo weakly in L\Rl : W). (3.9) 

Using p.7p and the two convergence properties p.8p and (|3.9p . we can proceed in the same 
way as Lemma 10 in |CJMTU] and prove that Uoo satisfies 



m — 1 



2 



dx ^ 0, (3.10) 



which imphes that Uoo is a radial stationary solution to (jl.l[) . hence is indeed in the family 



Next we will prove that u{-,t) Uoo uniformly in L°°(R'^) as t — > oo. In order to prove 
this, we make use of the monotonicity of the second moment of u{-, t) in time. By combining 
the following Virial identity 



d 
It 



■ / \x\^u{x, t)dx = 2{d - 2)F[u{t)] for all t (3.11) 
with the fact that the minimizcr of J- has free energy 0, it is shown in |BCL] that 



M2H-,t)] 



i |x|^w(2;, i)dt is non-decrasing in t. (3-12) 



This implies that any subsequence of w(-, i) can converge to only one limit: if not, then we can 
find some another sequence {t'^} increasing to infinity, such that uit'^) converges to another 
stationary solution u'^ uniformly as 7i — s- oo, where u'^ is also in the family (|1.4p . Since u[t[^) 
are uniformly bounded and uniformly compactly supported, we have M2[M(tJJ] — > M2[u'^]. 
On the other hand for the time sequence we have M2[u(i„)] — >■ M2[uoa\i hence ()3.12|) 
implies that and u'^ must have the same second moment. Since both u^o and u'^ are 
within the family ()1.4p . they can have the same second moment only if they are the same 
stationary solution. □ 

Remark 3.3. Since the proof is done by extracting a subsequence of time, we are unable 
to obtain the rate of the convergence. Moreover, it is unknown whether the limit Uoo has 
continuous dependence on the initial data uq. Wc also point out that the above proof is for 
radial solution only; for general initial data the difficulty lies in the fact that wc arc unable 
to bound the solution in some compact set uniform in time. 



3.2 Convergence towards self-similar solution for subcritical mass, 
radial case 

In this subsection, we prove that every radial solution with subcritical mass and compactly 
supported initial data would converge to some self-similar solution which is dissipating with 
the same scaling as the solution of the porous medium equation. 

Let u be the weak solution to (|1.1[) . with mass A e (0,A/c). Following ^ and [BCLj . we 
rescale u according to the scaling from porous medium equation: 

;i(A,r) = (t + l)u(x,t); A = a;(t + l)"^/''; t = ln(i + 1). (3.13) 
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Then /x(A,0) = u{x,0), and fi{X,T) solves the following rcscalcd equation 

fir = A/" + V • (A^vi^) + V ■ ifiVifi * M)). (3.14) 

It is pointed out in Theorem 5.2 of |BCLj that the free energy associated to the rescaled 
problem (\3.14\i is 

g{K;t)):^[ ('^^™ + i^(A^*^) + ^)dA, (3.15) 
jRd V TO - 1 2 2d J 

and for any mass A G (0, Mc), there is a unique minimizer fiA of G in Za subject to translation, 
where Za := {h G L^iW^) n L™(R'') : \\h\\i = Af and \x\'^h{x)dx < oo}. In addition, fiA 
is continuous, radially decreasing and has a compact support, and fiA satisfies 



+ i + l^iA:^ = (3.16) 



m — 1 dr ^ d aar''-^^ 
in its positive set, where the mass function M is as defined in p.2p . 

Since fiA is a stationary solution of p.l4p . if we go back to the original scaling, fiA gives a 
self-similar solution of (|1.1|) : 

7/^(x,t) = (t + l)-VA(^^^:^yY7d)' (3-17) 

It is then asked in |BCL| and |B] that whether this self-similar solution attracts all global 
solutions. 

We will first prove that all radial solutions to the rescaled equation p.l4p converge to fj,A- 
The following lemma construct a family of explicit solutions to p.l4p . which all converge to 
fjLA exponentially fast as t — >■ oo. 

Lemma 3.4 (A family of explicit solutions). Suppose d > 3 and m = 2 — 2/d. For 
< A < Mc, we denote by fiA the stationary solution of p.l4p . Let p, be defined as 

where R(t) solves the ODE 

^ ^ d^R'^ ' (3.19) 
i?(0) - i?o, 

where Rq > is a constant. Then for any Rq > 0, /2(A, r) is a weak solution to p.l4p . 

Proof. Since p, is a self-similar function, it can be easily verified that /2 solves the following 
transport equation 

p, + v.(m|<I>a).o, 

On the other hand, note that p. 141) can also be written as a transport equation 

//r = V • {fiv), 



12 



where 

a- ^ V7..m-i , ^ , i\/(|A|,r;/^) A 
+d+ adlAI'^-i |A|' 

Therefore, to prove that p. solves p.l4p . it suffices to verify that 

- ^ f M'"- + for < . < (3.20) 

Since /x is a rescaling of , 



where in the last inequality we used the fact that m is the critical power, i.e. m = 2 — 2/d. 
For T2 in p.20p . the definition of p. gives 

Now recall that satisfies p.l6p in its positive set, which implies 

RHSofS = -Vt ( ^#A^r^(i^) + + "^iT^ri 1 + ^(1 - l^)r 

Rjr) 

where the last equality comes from the definition of R in p.l9p . This verifies that p.20p is 
indeed true, which completes the proof. □ 



Next we use the family of explicit solution constructed above as barriers, and perform mass 
comparison between the real solution and the barriers. 

Proposition 3.5. Suppose d>3 and m = 2 — 2/d. Let fi{X,T) be a radial weak solution to 
(|3T4|) with mass < A< Mc, where the initial data 0) S L]_{R'^; (1 + \x\'^)dx) H L°°(M'') 
is continuous and compactly supported. Then as t ^ 00, the mass function of ^ converges to 
the mass function of exponentially, i.e. 

sup |M(r, r; /i) — M{r, r; ^a)\ < Ce~^ , 

r 

where fiA *s o,s defined in p.l6p . and C depends on d,A and /i(-,0). 



Proof. Without loss of generality we assume that /i(0, 0) > 0. (When /i(0, 0) — 0, from the 
same discussion in l3.1[ /i(0, r) will become positive after some finite time.) Then we can find 
Rqi sufficiently small and i?02 sufficiently large, such that 

^A..(^)^M-,o)^^/..(— ), 
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where in the first inequahty we used that /i(0, 0) > 0, and in the second inequality we used 
||At(-,0)||oo < oo. 

Let /ii(A,r) and /i2(A, r) be defined as in (|3.18p . with i?(0) equal to i?oi and R02 respec- 
tively. Then Lemma [3.41 says that both /ii and fj,2 are solutions to (|3.14p . Note that p.l4p 
is a special case of (|2.ip , hence the mass comparison result in Proposition 12.41 holds here as 
well, which gives 

yU2(-,T) -< h{-,t) -< ^ll{■,T) for all r > 0, 

or in other words, 

M(-,T;/i2) < M{-,t;h) < M{-,t;hi) for aU r > 0. 

It remains to show that 

svLp\M{r,T; - M{r; ha)\ < Ce"^ for i = 1,2. 

r 

Recall that both ^^'s are scalings of fiA with scaling coefficient Ri{r)^ hence 

\M{r,T■^i,)-M{r■^lA)\ - \M{-^- ^iA) - M{t- ^iA)\■ (3.21) 

Since ^a is bounded and compactly supported, it suffices to show that Ri{T) — )■ 1 exponen- 
tially as ?■ — > 00. Recall that Ri = ^(-^ — l)i?i with initial data i?oi for i = 1,2, a simple 

calculation reveals that \Ri{T) — 1| < CiC^'^, where Ci depends on Rqi. This implies that the 
right hand side of p.2ip decays like e~'^, which completes the proof. □ 

Making use of the explicit barriers /xi and fi2 constructed in the proof of Proposition l3.5| we 
get exponential convergence oi fj,/A towards the /x^/A in the p-Wasserstein metric, which is 
defined below. Note that the Wasserstein metric is natural for this problem, since as pointed 
out in |AGS| and |CMV| . the equation (|l.ip is a gradient flow of the free energy (|1.2p with 
respect to the 2- Wasserstein metric. 

Definition 3.6. Let jii and /i2 he two (Borel) probability measure on M'' with finite p-th 
moment. Then the p- Wasserstein distance between fii and ^2 is defined as 



p(/^i,M2):=( inf I / \x~y\'PT:(dxdy)\Y , 

where V{[1\t[12) is the set of all probability measures on x M.'^ with first marginal /xi and 
second marginal fi2- 

Corollary 3.7. Let d > 3, and to = 2 — |. Let /i(A,T) and jiA be as given in Proposition 
\3.5l Then for all p > 1, we have 



WA^,^)<Ce-^, 



where C depends on d and /i(-,0). 

Proof. See the proof of Corollary 5.8 in [KYj . □ 
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Rescaling back to the original space and time variables, we have 

Thus Corollarv 13 . 71 immediately yields the algebraic convergence towards the dissipating self- 
similar solution (|3.17p : 

Corollary 3.8. Let u(x,t) be a radial solution to (jl.ip with mass < A < Mc, where the 
initial data u{-, 0) G L^{M.'^] {l + \x\'^)dx)r\L°° {R'^) is continuous and compactly supported. Let 
UA be the dissipating self-similar solution with mass A defined in p.l7|) . Then u/A converges 
to ua/A in Wasserstein distance algebraically fast as t oo. More precisely, 

w,C^,^)<ct-^'~''i\ 

where C depends on d^A and u(-,0). 



3.3 Convergence towards self-similar solution for subcritical mass, 
non-radial case 

In this subsection, we consider the rescaled equation p.l4p with general (possibly non-radial) 
initial data. The key result here is that when the mass A < Mc is sufBcicntly small, the radially 
symmetric stationary solution jj,A as defined in (|3.16p is the unique compactly supported 
stationary solution (in rescaled variables). Then a similar argument as in Theorem 13.21 shows 
that every solution to (|3.14[) with small mass and compactly supported initial data converges 
to Ha- After scaling back to the original variables, we immediately obtain the convergence 
towards the self-similar solution if the mass is small. 

We first prove a L°°-regularization result, saying that if the initial mass is small, then the 
L°° norm of solution to p.l4p will become small after unit time, regardless of the L°° norm 
of the initial data. We point out that a similar L°°-regularization result is proved in [SSj for 
the 2D case with linear diffusion, using a De Giorgi type method. 

Lemma 3.9. Suppose d > 3 and rn = 2 — 2/d. Let /i(A, r) he a weak solution to (|3.14p with 
mass < j4 < Mc/2, where the initial data fig E L^(M'^; {l + \x\'^)dx)r\L°° (R'^) is continuous. 
Then we have 

||m(-,^)IIoo <KA:^CA^/''forallT>l, (3.22) 
where C is some constant depending only on d. 

Proof. Similar argument as the proof of Lemma 13.11 yields that 

/!*(•, r) ^ ^(^r) for all r > 0, (3.23) 

where /i(-, r) is the solution to ()3.14p with initial data /Iq. Since is radially symmetric and 
bounded above, we can find i?o sufficiently small, such that -< fiAii^) , '^hcre /i^ is as 
defined in (|3.16p . It then follows from Proposition 12.41 and Lemma 13.41 that 

^ 1^^'^(i?77y) " ^ 0' (3.24) 
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where R{t) satisfies the ODE ((XTOl) with initial data i?(0) Rq. Combining ((X^ and 
(13211), we obtain that 

||/^(-,t)||oo = ||ai*(-,t)||oo < -^^IImaIIoo for ah r > 0. 

In order to bound tlie riglit hand side of the above inequahty, we first find an upper bound 
for 1/R{t)'^. It can be readily verified that R{t) = min{iT~, ^} is a subsolution to p.l9p 
for any i?o > 0, which implies that R{t) > R{t) > i for all r > 1, thus -^^^p < S'^ for all 
r > 1. 

Next we will estimate ||/iA||oo- Note that fiA is radially decreasing for any < A < Mc, 
moreover ||/xa||oo = A'a(O) is increasing with respect to A. Therefore we readily obtain a 
rough bound ||^a||oo < Ci for all < A < Mc/2, where Ci — ^a/<:/2(0) only depends on d. 

Note that this rough boimd of ||/iA||oo gives us an upper bound for the velocity field given 
by the interaction term, namely 

dA,A*M) = ^^^<^. (3.25) 

To refine the bound for |l/iyi||oo: we compare fiA with fl^, where (Ia is the radial stationary 
solution to the following equation 

M. = A^" + V.(Mvii±§l^). (3.26) 

Making use of (|3.25p . mass comparison yields that fiA -< fJ-A, which implies /^a(0) < /iA(0). 
On the other hand note that (|3.26p is a Fokker-Planck equation, whose stationary solution is 
given by 

. / (l + Ci)(m-l) ,^„V/('»-i) 

^^=1^-^ 2d^ + 

where Ca > is the unique constant such that UmaHi = A. A simple algebraic manipulation 
shows that /iA(0) < CA^/'^, where C > depends only on d, therefore we can conclude. □ 



The next lemma shows that if the mass is small, any solution with compactly supported 
initial data will eventually be confined in some small disk. 

Lemma 3.10. Suppose d > 3 and m = 2 — 2/d. Then for any Rq > 0, there exists some 
sufficiently small Aq > 0, such that all weak solutions to p.l4p with continuous and compactly 
supported initial data and mass < A < Aq will be eventually confined in B{0,Ro). 



Proof. Let /^(A, r) be a weak solution to (|3.14l) with continuous and compactly supported 
initial data and mass < A < Aq, where Aq is a small constant depending on Rq and d to 
be determined later. 

In the proof of this lemma we take r = 1 to be the starting time, in order to take advantage 
of the estimate p.22[) . Our goal is to show that if the support of /^(•, 1) is contained in some 
disk B{0, R) where R > Rq — Ka and Ka is as defined in p.22p . then there exists some time 
T > 1 to be determined later, such that 

supp /^(•,t) C B{0,R + Ka) for all r G [1,T], (3.27) 
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moreover at time T the support can be fit into some disk smaller than B{0, R)^ namely 

supp n{-,T) C B{0,R- Ka/2). (3.28) 

By taking T as the starting time and repeating this procedure, we know that eventually the 
support will be confined in B(0,Rq). 

In order to deal with non-radial solution, we shall construct barriers in the density sense 
instead of in mass sense. Although comparison principle in density sense does not directly 
hold for p.l4|) due to the nonlocal term, if we treat V{X,t) := /i * TV as a fixed a priori 
potential, then (|3.14p becomes 

M. = A/i'" + V. (^A^V(^ + nA,r))), (3.29) 

which is a porous medium equation with a drift, and the weak solutions to it enjoy the 
comparison principle due to |BH] . It follows from p.22p that the following estimates of V 
holds: 

Ay(A, r) < sup ^ < Ka for a e M'', t > 1, 

and 

|VF(A,r)| < sup( fi* J < C(d)Al-^ for XeW^,T>l. 

Note that in both estimates above, the right hand side will go to zero as ^ — > 0. We also 

3 _ 2 I \ |2 

point out that ii R^ A'^ d^, then W will be dominated by around r = R. 

Next we will construct some explicit supersolution p. to p.29p . More precisely, we hope to 
find a continuous radially decreasing function /i defined in {r > R — Ka} x [1, T] for some T, 
such that fl satisfies the following inequality 

d - 1 ilr 3 2 

Mr > drrfl"' + {d,- + )( V) + + C(d) A rf " for all r > i? - Ka,t e [1, T], 

r d 

(3.30) 

while also satisfies the initial condition 

/i(r, 0) > Ka for aU R- Ka <r < R, (3.31) 
and the boundary condition 

fl{r, t)>Ka atr = R- Ka for all t e [1, T]. (3.32) 
The inequalities p.30P " P.33|) guarantees that fi is a supersolution to p.29|) . If A is small 

3 _ 2 

enough such that R > CA-^ for some large constant C depending on d, one can check that 

/1(A, r) = [2Ka - r{r - {R - Ka))] 

satisfies the inequalities (|3.30p ~ p.32p for 1 < r < 4, hence comparison principle yields that 
fl < (lin {r > R - Ka} for ah r £ [1, 4]. 

The reason we choose fl as above is that its support will shrink after some time: note that 
its support stays in B{0, R + Ka) for r G [1, 4], and most importantly, at r = 4, the support 
of fl can be fit into a disk smaller than B{0, R), namely 

supp /2(-,4) C B{0,R- Ka/2). (3.33) 
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Since comparison property gives that supp r) C supp fl{-, r) for all r G [1,4], we immedi- 
ately obtain p.27p and (|3.28p . which complete the proof. □ 



Making use of the above two lemmas, in the next theorem we show that when the mass is 
sufficiently small, there cannot be any non-radial stationary stationary solutions. 

Theorem 3.11. Suppose d > 3 and m = 2 — 2/d. Then when Q < A < Afc/2 is sufficiently 
small, the compactly supported stationary solution to (|3.14p is unique. 



Proof. Due to Corollarv 13.71 we know that for any < A < Mc, there does not exist any 
compactly supported radial stationary solution other than ^a- Hence it suffices to prove 
that when A is sufficiently small, every compactly supported stationary solution is radially 
symmetric. 

Suppose VAi^) is a compactly supported stationary solution to (|3.14p . which is not radially 
symmetric. Since va is stationary, it satisfies 



TO — 1 a 



VA*J^+^^C in {vA> 0}, (3.34) 



where different positive components of va may have different C's. Heuristically, the idea is 
to argue that the term va* -N" must be more "roundish" than i^^v^' if VA is non-radial, 
thus get a contradiction. 

We point out that p.34p implies that va is continuous in and smooth inside its positive 
set. This enables us to find two points a, 6 e M'* in the same connected component of {va > 0}, 
satisfying \a\ = \b\ and 

VA{a) - VA{b) = sup {va{x) - VA{y)) > 0. (3.35) 

|a:| = |y 

We claim that when A is sufficiently small, the following inequality holds 

-^y^-\a) ~ v'r\b)\ > \{vA*m'^)~i'^A*mb)l (3.36) 

TO — 1 ' ' 

then (|3.36p would contradict p.34[) . 

We start with the left hand side of p.36|) : Lemma [3.91 implies that both VA{a) and VA{b) 
are much smaller than 1 when A is small. Since 0<to — 1<1, it follows that 



TO 



j\v^-^ia)-v^r^{b)\>\vA{a)-VAib)\. 



TO 

if A is sufficiently small. In order to prove p.36p . it suffices to show that 

\ivA*M)ia) ~ {vA*M){b)\ < \vA{a) - VA{b)\. (3.37) 



We introduce a linear transformation T : R'^ — > R'^ which is a rotation that maps a to b. 
Then radial symmetry of Af yields that {va * A/") (6) = {{va oT) * Af){a). 
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In addition, T being a rotation implies that |r(a;)| = |a;| for any x = W^, hence from the 
way we choose a and 6, we have \iyA{T{x)) — va{x)\ < i^Aia) — VA{b) for any x S IR''. Thus 

\{vA*J^){a) - ii^A*J^){b)\ = \{iyA*J^){a) - iii^A°T)*N'){a)\ 

< [ \i^Aiy)~iyA{T{y))\\Mia~y)\dy 

< {vA{a)-VA{h)) [ \M{y)\dy, 

Jb{o,r) 

where B{0,R) is the smaUest disk that contains the support of i^a- Now we make use of 
Lemma I3.10[ which shows that we can fit the support of va into an arbitrarily small disk by 
letting A be sufficiently small. Therefore we can choose R such that J^^^ \N'{y)\dy < 1/2, 
then let A be sufficiently small such that supp i/a C B{0,R). This gives us p.37|) . which 
leads to a contradiction and hence completes the proof. □ 

Remark 3.12. For general < ^ < Mc, we arc unable to prove the uniqueness of the 
compactly supported stationary solution. The difficulty lies in the fact that for larger mass 
we are only able to show the support lies in a disk with radius 0(1). Hence instead of p.37p . 
we can only obtain {{i/a *N'){a) — {i^a *Af){b)\ < C|;^A(a) — f^Aib)], where C might be a large 
constant, which stops us from getting a contradiction. 

Once we obtain the uniqueness of compactly supported stationary solution for small mass, 
the following corollary shows that all solution with compactly supported initial data must 
converge to this unique stationary solution as r — > cx). 

Corollary 3.13. Suppose d > 3 and m ~ 2 — 2/d. Let /i(A, r) be a weak solution to p.l4p 
with mass < A < Mc/2 being sufficiently small, where the initial data IJ.{-, 0) G Li'|_(M'*; (1 + 
|a;p)(ia;) H L°°{W^) is continuous and compactly supported. Then as t oo, we have 

M■,T)~^,A{■)\\oo~^0, (3.38) 

where fiA is as defined in (|3.16p . 

Proof. The proof is similar as the proof of Theorem 13.21 and actually it is simpler here since 
there is a unique stationary solution, instead of a family of stationary solution in the case of 
Theorem [321 

When the initial data 0) is bounded and compactly supported. Lemma 13.91 and Lemma 
13. 101 shows that t) would be uniformly bounded and stay in some fixed compact set for all 
r > 1. In addition, the continuity result in [D] indicates that //(A,t) is uniformly continuous 
in space and time in R'^ x [1, oo). 

As a result, for any time sequence r„ that increases to infinity, using the same argument as 
in the proof of Theorem 13. 2| we can extract a subsequence r„^ such that fi{-,T„^) uniformly 
converges to some continuous function Hoo, where fi^o is a compactly supported stationary 
solution. Theorem 13.111 ensures that /Xoo niust coincide with (ia when A is sufficiently small, 
yielding that r) indeed converges to /i^ uniformly as t — > oo. □ 
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Remark 3.14. Since fi{-,0) is confined in some compact set for all time, (j3.38p implies that 
\\fi{-,T) — fiA{-)\\p as T ^ oo for all p > 1. Now if we scale back to the original variables, 
it immediately follows that — w^(-)||p — )■ as < — 5- cx) for all p > 1, where ua is the 

dissipating self-similar solution as defined in p.l7p . However the rate of convergence here is 
unknown, since the proof is done by extracting a subsequence of time. 

4 Application to aggregation models with repulsive-attractive 
interactions 

In this section we consider the foUowing integro-diffcrcntial equation 

ut = V ■ (uVK *u), (4.1) 

where the interaction kernel K has a repulsion component in the form of the Newtonian po- 
tential M{x) = — (t;_2)er|j|a;|<J-ii an attraction component satisfying the power law, namely 

K{x)=M{x) + ^\x\'^, (4.2) 

where 2 — d < q < 2, and when q = the second term is replaced by In |a;|. 

The global existence of weak solution is established in |FH| for q > 2 — d. Next we show 
that mass comparison holds for (|4.ip between weak solutions. 

Proposition 4.1. Let ui{x,t), U2{x,t) be two radially symmetric weak solution to ()4.ip . 
which are compactly supported for all i > 0. //ui(-,0) < W2(-,0), then ui{-,t) -< U2{-,t) for 
all t > 0. 

Proof. Note that A^|a;|'' = {q + d — 2)|a;|'~^, which is locally integrable in M''. nonnegative 
and radially decreasing when 2 — d < q < 2. Therefore the conditions (C), (Kl), (K2'), 
(VI) are met, (where (K2') is as defined in Remark l2.5p . and (|4.ip becomes a special case 
of the general equation (|2.ip in Section 2. Due to the discussion in Remark l2.51 we can apply 
mass comparison to the compactly supported solutions. 

Without loss of generality we assume that both wi(-, 0) and U2(-, 0) are continuous, and for 
general initial data we can use approximation. Due to Theorem 2.5 of [FH] , we have Ali is 
in both space and time for all t, where i = 1,2. Now we can apply Proposition 12.41 to 
conclude that wi(-, t) -< U2{-,t) for all time: although the proof of Proposition 12 .41 requires Mi 
be in space, the requirement are only for the diffusion term. Since the right hand side 
of (|4.ip only has aggregation terms, continuity of Mi is sufficient. □ 

The following existence and uniqueness result of a stationary solution is established in }FH) . 

Proposition 4.2 ( |FH| . Theorem 3.1). For every q > 2 — d and mass A > 0, there exists 
a unique radius Ra (that depends on q and d only) and a unique steady state Ug of the 
aggregation model (|4.ip - (|4.2p that is supported on B{0,Ra), has mass A and is continuous 
on its support. 
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In addition, Ug is radially decreasing in B{0,Ra) if q < 2, and is radially increasing if 
q>2. 



It is proved by [FHj that in its positive set, Ug satisfies VK * Ug = 0, which can also be 
written as —Us + A{^\x\'^) * = 0, i.e. 

- Us + {q + d ~ 2)\x\'^-^ * Us ^ in {u, > 0}. (4.3) 

In the proposition below, we construct a family of explicit subsolutions, all of which are 
compactly supported and converge to the stationary solution Us exponentially fast. 

Lemma 4.3 (A family of explicit subsolutions). Suppose d > 3 and 2 — d < q < 2. Let 

Us be the stationary solution to (j4.ip with mass A, as given by Proposition \4-S\ We define 
the self-similar function u as 

^-.(^), (4.4) 

where R{t) solves the ODE 

i R{t) ^ Cil - R'+^~')R-''+' 
\ R{0) = Ro, 

where Rq > 1, and Ci is some fixed constant only depending on q,d and A. Then for all 
t G [0, oo), u{x,t) is a subsolution to (|4.ip in the mass comparison sense. 



Proof. The proof is similar to the proof of Lemma 13.41 Since w is a self-similar function, it 
can be easily verified that u solves the following transport equation 

^* + V.(f.^.)^0, 

which implies that the mass function of u satisfies 

Mt{r,t;u)=adr''''Mrir,t;u) (-|^)r. (4.6) 

Our goal is to show that tl is a subsolution to ()4.1|) in mass comparison sense, which is 
equivalent to the following inequality due to Lemma 12.31 



Mt{r,t;u) < adr'^-^Mr{r,t;u) 



M{r,t;u) ^ M{r,t;u) 



(4.7) 



where M(r, i; It) = Jg^g ^^(g + d - 2)(|a;|* * u){x,t)dx. 



We point out that is nonnegative by definition, since Mr{r,t;u) ~ adr'^u, where u is 
nonncgativc. By comparing (j4.6|) and (j4.7p . it suffices to prove that the following inequality 
holds for u(r, t) for ah < r < R{t): 

m ^ M{r,t-u) M{r,t-u) 

R{t) - OdT'^-^ ^ ' ^^■''^ 
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Next we will investigate the terms on the right hand side of (|4.8p . Since u is defined as a 
continuous scaling of Ug, for the first term on the right hand side of (|4.8p . we have 

For the second term, we obtain that for any < r < R{t), 

M{r,t;u) = iq + d-2) [ [ \x ~ y\^-^ -^Usi-^)dydx 

JB(0.r) JR-i R{t) R{t) 

= {q + d-2)R{tY-^ I I \^ z\i-^u,[z)dzdx {z:=y/R{t)) 



iq + d-2)R{tr-' f i\xr^,u,){-^)dx 

mr-^ I ^s{^Mx (by m) 

JB(Q,r) ^{'') 



'-8(0, r) 

Putting the above two equations together yields 

RHS of g^) = (-1 + R{tf+i-^) ^ 

> Ci{-1 + R{tf+''-^)R{t)-'^r. (4.9) 

Where Ci only depends on q, d and A. Here in the last inequality we used the fact that Ug is 
radially decreasing, which implies that M{r;Us) > CA\B{0,r)\ for all < r < Ra, where Ca 
is the average density of Ug in its support -6(0, Ra)- 

Since R = Ci{l — by definition, the above inequality implies that ()4.8p is 

true, which completes the proof. □ 

Remark 4.4. Similarly, we can construct a family of explicit supers olutions u in the mass 
comparison sense: here u is defined in ()4.4p . where R{t) solves the ODE ()4.5p with initial 
data < Rq < 1, and the constant Ci in (j4.5p is replaced by some other fixed constant C2, 
which also only depends on q, d and A. 



Making use of the subsolutions and supersolutions we constructed in Lemma 14.31 and Re- 
mark |4]4] respectively, we next prove that all radial solutions with compactly supported initial 
data will converge to the unique stationary solution exponentially fast. 

Theorem 4.5. Suppose d > 3 and 2 — d < q < 2. Let u be a weak solution to (j4.ip with 
initial data uq and mass A, where uq £ L"'^(IR^) H -L°° (M'^) is non-negative, radially symmetric 
and compactly supported. In addition, we assume that uq is strictly positive in a neighborhood 
0/0. Let Ug be the unique stationary solution with mass A, as given by Proposition \4-2\ Then 
as t — > 00, u{-,t) converges to Ug exponentially fast in Wasserstein distance. 



Proof. The proof is similar to the proof of Proposition l3.5l Since uq is bounded, and is strictly 
positive in a neighborhood of 0, we can find Rqi sufficiently large and -R02 sufficiently small. 
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such that 

1 • , 1 • , 

Let Ml be the subsolution given by Lemma 14.31 with initial data -j^Usl-n—, 0), and U2 be the 
supersolution given by Remark l4.4l with initial data -^Us{-n—,0)- Then mass comparison in 
Proposition 12.41 yields ui{-,t) -< u{-,t) -< U2{-,t) for all t. 

Now it suffices to show that ui and U2 both converges to Ms exponentially fast in Wasserstein 
distance. Note that ui{-,t) is a continuous scaling of m^ with scaling coefficient R{t), where 
R{t) satisfies the ODE (|4.5p and hence converges to 1 exponentially. More precisely, we have 



\R{t) - 1| < C[e 



-Ci{d+q~2)t 



where Ci is as given in Lemma 14.31 and C[ depends on q^d and i?oi- Similar result hold for 
the supersolution U2- Then argue as in Corollarv l3.71 we have for all p > 1 that 

W,C^,^)<c,e-^^^'+^-^^' for ^^ 1,2, 

where Ci depends on i?oz respectively. Since m^ is squeezed between ui and U2 in the mass 
comparison sense, the inequality above yields 

Wp{^^, ^) < C'e-^(''+«-2)* for ^ = l,2, 

where C :~ min(Ci,C2) depends on d,q and A, while C" depends on d,q,A and m(-,0). □ 



References 

[AGS] L. Ambrosio, N. Gigli and G. Savarc, Gradient flows in metric spaces and in the 
space of probability measures. 2nd ed. Lectures in Mathematics, ETH Ziirich. Basel: 
Birkhaiiser., 2008. 

[BCLR] D. Balague, J. A. Carrillo, T. Laurent and G. Raoul, Nonlocal interactions by 
repulsive-attractive potentials: radial ins/stability, preprint, arXiv: 1109.5258, 2011. 

[BK] J. Bedrossian and I. Kim, Global existence and finite time blow-up for critical Patlak- 
Kcllcr-Segel models with inhomogeneous diffusion, preprint, arXiv:1108.5301, 2011. 

[BRB] J. Bedrossian, N. Rodriguez and A.L. Bertozzi, Local and global well-posedness for 
aggregation equations and Patlak-Keller-Segel models with degenerate diffusion, sub- 
mitted. 

[BH] M. Bertsch and D. Hilhorst, A density dependent diffusion equation in population 
dynamics: stabilization to equilibrium, SIAM J. Math. Anal., 17(1986): 863-883. 

[BKLN] P. Biler, G. Karch, P. Laurengot and T. Nadzieja, The 87r-problem for radially sym- 
metric solutions of a chemotaxis model in the plane, M3AS, 29 (2006): 1563-1583. 

[B] A. Blanchet, On the parabolic-elliptic patlak-keller-segel system in dimension 2 and 
higher, preprint, arXiv:1109.1543, 2011. 



23 



[BCL] A. Blanchct, J. Carrillo and P. Laurcngot. Critical mass for a Patlak-KcUer-Segel 
model with degenerate diffusion in higher dimensions, Calc. Var., 35(2009): 133-168. 

[BCM] A. Blanchet, J. A. Carrillo and N. Masmoudi, Infinite time aggregation for the critical 
Patlak- Keller-Segel model in Comm. Pure Appl. Math., 61(2008): 1449-1481. 

[BDP] A. Blanchet, J. Dolbeault and B. Perthame, Two-dimensional Keller-Segel model: 
Optimal critical mass and qualitative properties of the solutions, E. J. Dijf. Eqn, 
44(2006): 1-32. 

[BEDF] A. Blanchet, J. Dolbeault, M. Escobedo and J. Fernandez, Asymptotic behavior 
for small mass in the two-dimensional parabolic-elliptic Keller-Segel model, J. Math. 
Anal. Appl, 361(2010): 533-542. 

[BDF] M. Burger, M. Di Francesco and M. Franek. Stationary states of quadratic diffusion 
equations with long-range attraction, i arXiv.l 1 03. 5365\ 2011. 

[CJMTU] J. A. Carrillo, A. Jiingel, P.A. Markowich, G. Toscani and A. Unterreiter, En- 
tropy dissipation methods for degenerate parabolic problems and generalized Sobolev 
inequalities, Monatsh. Math., 133 (2001): 1-82. 

[CMV] J. A. Carrillo, R. J. McCann and C. Villani, Contractions in the 2-Wasserstein length 
space and thermalization of granular media. Arch. Rat. Mech. Anal. 179(2006): 217- 
263. 

[D] E. DiBenedetto, Continuity of weak solutions to a general porous media equation, 
Indiana Math. Univ. J., 32 (1983): 83-118. 

[DP] J. Dolbeault and B. Perthame, Optimal critical mass in the two dimensional Keller- 
Segel model in R^, C.R. Acad. Set. Paris, Ser I Math, 339(2004): 611-616. 

[FH] R. C. Fetecau and Y. Huang, Equilibria of biological aggregations with nonlocal 
repulsive-attractive interactions, preprint, arXiv:1109.2864, 2011. 

[FHK] R. C. Fetecau, Y. Huang and T. Kolokolnikov, Swarm dynamics and equilibria for a 
nonlocal aggregation model, Nonlinearity, 24(2011): 2681-2716. 

[GGL] I. Gallagher, Th. Gallay and P.-L. Lions, On the uniqueness of the solution of the 
two- dimensional Navier-Stokes equation with a Dirac mass as initial vorticity. Math. 
Nachr., 278(2005): 1665-1672. 

[JL] W. Jager and S. Luckhaus. On explosions of solutions to a system of partial differntial 
equations modelling chemotaxis. Trans. Amer. Math. Soc, 329(1992): 819-824. 

[KS] E. F. Keller and L.A. Segel, Model for chemotaxis. J. Theor. Biol., 30(1971): 225-234. 

[KY] I. Kim and Y. Yao, The Patlak-Keller-Segel model and its variations: properties of 
solutions via maximum principle, to appear in SIAM J. Math. Anal, arXiv:1102.0092, 
2011. 

[KSUB] T. Kolokolnikov, H. Sun, D. Uminsky and A. L. Bertozzi, A theory of complex 
patterns arising from 2D particle interactions, Phys. Rev. E, Rapid Communications 
84(2011) 015203(R). 



24 



[H] D. Horstmann, From 1970 until present: the Keller-Segel model in chemotaxis and its 
consequences. /, Jahresber. Deutsch. Math.-Verein., 105(2003): 103-165. 

[LTB] A. J. Leverentz, C. M. Topaz and Andrew J. Bernoff, Asymptotic dynamics of 
attractive-repulsive swarms, SIAM J. Appl. Dyn. Syst., 8(2009): 880-908. 

[ME] A. Mogilner and L. Edclstein-Kcshct. A non-local model for a swarm. J. Math. Biol., 
38(1999): 534-570. 

[P] C. S. Patlak. Random walk with persistence and external bias. Bull. Math. Biophys., 
15(1953): 311-338. 

[SS] T. Senba and T. Suzuki, Weak solutions to a parabolic-elliptic system of chemotaxis, 
J. Func. Analysis, 47 (2001), 17-51. 

[TB] C. M. Topaz and A. L. Bertozzi, Swarming patterns in a two-dimensional kinematic 
model for biological groups, SIAM J. Appl. Math., 65(2004): 152-174. 

[TBL] C. M. Topaz, A. L. Bertozzi and M. A. Lewis, A nonlocal continuum model for bio- 
logical aggregation. Bull. Math. Biol, 68(2006): 1601-1623. 

[V] J. Vazquez, The Porous Medium Equation: Mathematical Theory, Oxford University 
Press, 2007. 



25 



